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Abstract 

We estimate the probability that a hard nucleon-nucleon collision is able to nucleate 
a seed of quark-gluon plasma in the surrounding hot and dense hadronic matter 
formed during a central collision of two large nuclei at AGS energies. The probability 
of producing at least one such seed is on the order of 1-100%. We investigate the 
influence of quark-gluon plasma formation on the observed multiplicity distribution 
and find that it may lead to noticeable structure in the form of a bump or shoulder. 
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1 Introduction 



One of the mysteries of heavy ion physics at Brookhaven National Laboratory's 
AGS is: // hadronic cascade event simulators like RQMD /fj/ and ARC /||/ produce 
energy densities approaching 2 GeV/fm 3 , yet agree with experiment, where is the 
quark-gluon plasma? After all, numerous estimates of the onset of quark-gluon 
plasma agree that it should occur at about that energy density, and if there is a first 
order phase transition, then the onset of the mixed phase would occur at an even 
lower density. One possibility is that no phase transition occurs even at these high 
densities, but it is difficult to understand how composite objects like hadrons can 
overlap so strongly in position space without the matter undergoing some qualitative 
change in character. A second possibility is that the distribution of observed hadrons 
in the final state is insensitive to the dynamics of the matter when it is most hot and 
dense. (Unfortunately there are no measurements of direct photons or dileptons at 
the AGS which might probe this stage of the collision.) There is some evidence for 
this which comes from artificially modifying hadronic cross sections at high density 
[[|. It may be understood by recognizing that once a system reaches local thermal 
equilibrium it is basically irrelevant how it got there. 

Recently we proposed a third possibility ||: Most collisions at AGS energies 
produce superheated hadronic matter and are describable with hadronic cascade simu- 
lators, but in rare events a droplet of quark-gluon plasma is nucleated which converts 
most of the matter to plasma. We estimated the probability of this to occur, using 
homogeneous nucleation theory, to be on the order of once every 100 to 1000 central 
collisions of large nuclei. Our estimate was based on the probability that thermal 
fluctuations in a homogeneous superheated hadronic gas would produce a plasma 
droplet, and that this droplet was large enough to overcome its surface free energy 
to grow. In this paper we consider another source of plasma droplet production 
which is essentially one of nonthermal origin. Specifically, we estimate the proba- 
bility that a collision occurs between two highly energetic incoming nucleons, one 
from the projectile and one from the target, that this collision would have produced 
many pions if it had occurred in vacuum, but because it occurs in the hot and dense 
medium its collision products are quark and gluon fields which make a small droplet 
of plasma. Although there is a large uncertainty in our estimates, we find that 
this inhomogeneous nucleation of plasma may be more probable than homogeneous 
nucleation by one to two orders of magnitude. 

In this paper we also consider the problem of observation of the effects of nu- 
cleation of plasma in rare events. We are guided by observations of multiplicity 
distributions in pp collisions at the CERN and Fermilab colliders. In those distribu- 
tion, one sees a shoulder developing at high multiplicity at an energy of 540 GeV, 
which turns into a noticeable bump at higher energies. The real cause of this struc- 
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ture is not known, but may be due to minijet production. If plasma is nucleated in 
some fraction of central nucleus-nucleus collisions at the AGS, a similar structure 
may develop. 

2 Kinetic Model of Hard Nucleon-Nucleon Colli- 
sions 

In this section we develop a simple kinetic model which allows us to estimate the 
number of high energy nucleon-nucleon scatterings occurring in the high density 
medium formed during a collision between heavy nuclei. These scatterings occur 
when a projectile nucleon penetrates the hot and dense matter to collide with a 
target nucleon which has also penetrated the hot and dense matter. The energy loss 
of the colliding nucleons must be taken into account to obtain a reasonable estimate 
of the energy available for meson production in the nucleon-nucleon collision. 

To first approximation we can visualize the initial stage of a heavy ion colli- 
sion at the AGS in the nucleus-nucleus center-of- momentum frame as two colliding 
Lorentz contracted disks. See Figure 1. At time t = they touch; subsequently 
they interpenetrate, forming hot and dense matter in the region of overlap. During 
this stage, additional matter streams into the hot zone even as this zone is expand- 
ing along the beam axis. The nucleons streaming in undergo scatterings with the 
hot matter already present, degrading their longitudinal momentum and producing 
baryonic isobars and/or mesons. Finally, at time to = L/2vj, all the cold nuclear 
matter has streamed into the region of overlap, and expansion and cooling begins. 
Here, L is the nuclear thickness, v is the velocity in the center-of-momentum frame, 
and 7 is the associated Lorentz contraction factor. This is a very simplified picture 
of the early stage of the collision, but it seems to semi- quantitatively represent the 
outcome of both the ARC and RQMD simulations |TJ, [|]. 

We are interested in the possibility that an incoming projectile nucleon suffers 
little or no energy loss during its passage to the longitudinal point z inside the hot 
and dense zone where it encounters a target nucleon which also has suffered little 
or no energy loss. The energy available in the ensuing nucleon-nucleon collision, 
y/s, can go into meson production. Suppose that a large number of pions would 
be produced if the collision had happened in free space. Clearly, the outgoing 
quark and gluon fields cannot be represented as asymptotic pion and nucleon states 
immediately. The fields must expand and become dilute enough to be called real 
hadrons. If this collision occurs in a high energy density medium, the outgoing 
quark- gluon fields will encounter other hadrons before they can hadronize. It is 
reasonable to suppose that this "star burst" will actually be a seed for quark-gluon 
plasma formation if the surrounding matter is superheated hadronic matter. We 
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need a semi- quantitative model of this physics. 

A fundamental result from kinetic theory is that the number of scattering pro- 
cesses of the type 1 + 2 — > X is given by 

dt J d 3 x J ^y§/i(x,Pi,0 J f 2 (x,p 2 ,t) v 12 a 1+2 ^ x (s 12 ) . (1) 

Here v n is a relative velocity, 



v (Pi ' P2) 2 - 

U12 = ^ , (2) 

where pj denotes the four-momentum of nucleon « and E t = \Jpf + rn? N its energy. 
The fi are phase space densities normalized such that the total number of nucleons 
of type i is 

A differential distribution in the variable Y is obtained by replacing a with da/dY. 

For our purpose it is reasonable to represent the colliding nuclei as cylinders 
with radius R and thickness L. All the action is along the beam axis. We assume 
that the phase space distributions are independent of transverse coordinates x and y 
and of transverse momentum. Integrating over the cross sectional area of the nuclei, 
and counting only those collisions that occur within the hot zone, yields 

A^i+2^x = ttR 2 dt J dz J dP **^ 2 2z fi(z,p lz , t) f 2 (z,p2 Z , t) v 12 er 1+2 ^x(si2) ■ 

(4) 

Here there is a change in notation: fi(z,pi Z , t)/2n is the probability per unit volume 
to find a nucleon i with longitudinal momentum p iz at longitudinal position z at 
time t. The integration limits on z ensure that the collisions under consideration 
really occur in the hot zone; see Figure 1. The integration limits on t mean that we 
only count those collisions which occur before the system begins its cooling stage. 
The depth in the hot zone to which nucleon 1 has penetrated is di = (vt + z)/2, 
and the depth to which nucleon 2 has penetrated is d 2 = (vt — z)/2. We neglect 
the decrease in velocity of the nucleons as they travel through the hot zone. This 
is an acceptable approximation because in the end we are interested only in those 
nucleons which suffer a small energy loss in traversing the hot matter. 
We construct the phase space distribution as follows: 

H(x, N) = probability that the nucleon has momentum fraction x after 
making N collisions; 
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S(N, d) = probability that the nucleon has made N collisions after pen- 
etrating to a depth d; 

J2N=oH(x,N)S(N,d) = probability that the nucleon has momentum 
fraction x after penetrating to a depth d. 

The distribution functions are normalized to unity. 

f-H{x,N) = 1 (5) 
Jo x 

oo 

E S(N, d) = 1 (6) 

N=0 

The phase space density of nucleon % is then taken to be 
dr> ■ dr 00 

J^f l ( z ^t)= 1 r ki ^ £ H( Xl ,N t )S(N t ,d t ), (7) 

Zn X * Ni=0 

where Uq is the average baryon density in a nucleus, about 0.145 nucleons/fm 3 . As 
a check, we can compute the number of nucleons which have entered the hot zone 
as a function of time. 

NT*(t) = J ^±f(z,pu, t) e(dt) = 2nR^n vt (8) 

The step function fixes the limits on the z integration. The number of participating 
nucleons grows linearly with time, and at time to we get Nf 31 (to) = irR 2 Lno, which 
is the total number of nucleons in the nucleus. 

The number of elementary nucleon-nucleon collisions can now be expressed as 

_, 2 2 2 f to , f vt , f 1 dxi r 1 dx 2 , \ 

A/i+2-.x = 7 n dt dz / v 12 cr 1+2 ^x{si2> 

JO J-vt Jo X\ Jo x 2 

oo oo 

£ £ H(x u N ± ) S(N U d ± ) H(x 2 , N 2 ) S(N 2 , d 2 ) . (9) 

JVi=0 JV 2 =0 

Since the nucleons' velocities are antiparallel the velocity factor is 

XlPo x 2 p 

Vl2 = 1 = + 1 = , (10) 

^x\pl + m 2 N ^x 2 2 pl + m 2 N 

where po is the beam momentum in the center-of-momentum frame. 

The survival function S(N, d) is characterized by the mean free path A of nucle- 
ons in the hot and dense hadronic matter. For a dilute gas the inverse of the mean 
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free path is the sum of products of the cross section of the nucleon with the density 
of objects it can collide with. 

A" 1 = Wi (11) 

i 

Average particle densities, including baryons and mesons, were computed in ref. 
[|] for the hot and dense matter under consideration. A plot of the density as a 
function of beam energy is shown in Figure 2. Assuming an average hadron- nucleon 
cross section of 25 mb, we find A = 0.4 fm at a laboratory beam energy of 11.6 
GeV/nucleon. This is very short, and just emphasizes the physics we discussed in 
the introduction concerning hadronic matter versus quark-gluon plasma. 

We assume that the collisions suffered by the nucleons are independent and can 
be characterized by a Poisson distribution. 

s <^M®" exp H) (i2) 

Here d/X is the average number of scatterings in a distance d. 

The invariant distribution function H(x, N) describes the momentum degrada- 
tion of a nucleon propagating through the hot zone. This distribution function was 
introduced in the evolution model of Hwa ||. In this model the nucleon propagates 
on a straight line trajectory and interacts with target particles contained within a 
tube with area given by the elementary nucleon-nucleon cross section <Jnn- Cser- 
nai and Kapusta |J solved the resulting evolution equations and found that the 
invariant distribution function in this model is given by 

N ( N \ (— In tY^ 1 

H(x, N) = *Xr n ) - w)N ~ n ( n - 1)1 + (1 " W)N 6{X ~ 1} ' (13) 

The 5-function represents elastic and soft inelastic contributions to the evolution of 
the nucleon through the matter. The probability w is the ratio of inelastic to total 
nucleon- nucleon cross section. It corresponds to the probability that the nucleon 
scatters inelastically and therefore drops out of the evolution described by H; it 
is approximately 0.8 in free space. Csernai and Kapusta found that it reduces to 
about 0.5 for nucleons propagating through a nucleus. This value allowed them to 
obtain a good representation of data with beam energies in the range of 6-405 GeV. 
In our case the nucleon is propagating through hot and dense hadronic matter. We 
keep w as a free parameter since we don't know how the value of w changes due to 
the thermal excitations and the increased density. 

We are interested in the number of pion-producing nucleon-nucleon collisions 
with a relatively high center-of-momentum energy squared s. Our basic result from 
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this section is 

d^n ard ^222 , n f to , f vt , f 1 dx! f 1 dx 2 



ttR 7 n ai n (s) at dz / v\ 2 ois — s 12 ) 

as JO J— ui JO Xi JO £ 2 



oo oo 



J2 £ # K Ni) S(N ± , d x (z, t)) H(x 2 , N 2 ) S(N 2 , d 2 {z, t)) . (14) 

JVi=0 ^2=0 



Here o"i n is the inelastic nucleon-nucleon cross section, and y/su is the total energy 
in the nucleon-nucleon collision where the nucleons have momentum fractions X\ 
and %2- 



3 Meson Production Cross Sections 

A phase transition to quark-gluon plasma will become thermodynamically favor- 
able if the energy density is large enough. The corresponding phase boundary in 
the temperature/chemical potential plane was explored in ||. Until now we have 
only selected nucleon-nucleon scatterings in which the total available energy yfs is 
large. In addition, we need to specify what fraction of this energy goes into meson 
production. In this section we estimate the pion number distribution function P n (s), 
which is the probability of producing n pions in a nucleon-nucleon collision in free 
space. The pion number distribution function is linked to the cross section a n for 
producing n pions by 

P n (s) = a n (s)/cr in (s) . (15) 

Given P n (s) we can estimate the number of nucleon-nucleon collisions that would 
lead to the production of n pions as 

r AE o <i/V. hard 
N n = dsP n (s)^^. (16) 

J s min U>b 

The lower limit of integration is fixed by kinematics and the upper limit is deter- 
mined by the beam energy. 

We shall approximate the pion number distribution function P n (s) with a bino- 
mial J?| and choose the parameters of this binomial such that we have some rough 
agreement with experiment ||. 



P n {s) = [ max £ n (l-0 nmaJE ~ " (17) 




The maximum number of pions produced in a nucleon-nucleon collision is determined 
by kinematics. 



n. 



Integer (^~ 2m -)\ (i 8 ) 
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The parameter £ is related to the mean multiplicity by 



= ^ L = ^r~ (iW + + z< n ->) ■ ( 19 ) 

''■max ''■max * ^ / 

Here (n) is the average pion multiplicity averaged over pp, pn and nn collisions 
while (nip), (n~) and (n~ n ) represent the average negative pion multiplicity in 
those collisions. All average multiplicities are functions of s, of course. The factor 
of 3 is due to isospin averaging. 

Experimental data were compiled and parametrized in |J as 

(n~) = -0.41 + 0.79F(s) 
(n- n ) = -0.14 + 0.81F(a) 

(n~ n ) = +0.35 + 0.77F( S ). (20) 



The function F was introduced by Fermi |10 



F(s) ^-^f'\ (21) 

with s measured in GeV 2 . The parametrizations in (|20"D describe the data rather well 
except in the threshold region. We approximate the inelastic nucleon-nucleon cross 
section a m by the inelastic proton-proton cross section. A convenient parametriza- 
tion is given in |§, 

a in = 30.9 - 28.9pl 2A6 - 0.835 lnp L + 0.192 \n 2 p L , (22) 

where pi is the laboratory momentum in GeV/c and the cross section is in mb. This 
parametrization is good for pi > 0.968 GeV/c. 

The pion production cross sections, as described above, are displayed in Figure 
3. They have the right shapes and the right orders of magnitude compared to data 
[§]. However, direct comparison is not possible. First of all, data generally does 
not exist for final states with 7r + , 7r~, and 7r°. Usually, exclusive experiments can 
only measure charged mesons or neutral mesons, not both. Secondly, we have not 
been so sophisticated as to include vector mesons, the r\ meson, and kaons. For 
our purpose such sophistication is probably not necessary. We care only about the 
probability that a nucleon-nucleon collision leads to a significant amount of energy 
release in the sense of conversion of initial kinetic energy to meson mass. We are 
essentially basing our results on the total inelastic cross section, the average meson 
multiplicity, kinematics, and entropy. Our analysis would be better if we had a 
handle on the width of the multiplicity distribution, averaged over the initial state 
isospin and summed over the final state isospin. 
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4 Star Burst Probabilities 



In this section we put together the ingredients developed in the last two and compute 
the number of star bursts which may become nucleation sites or seeds for plasma 
formation and growth. 

The nucleon-nucleon collisions may be referred to as primary-primary, primary- 
secondary, and secondary-secondary, depending on whether the nucleons have scat- 
tered from thermalized particles in the hot zone (secondary) or not (primary). The 
easiest contribution to obtain is the primary-primary. All integrations and summa- 
tions can be done analytically with the result 



dN „™ _ {Xino 



ds \ w 



L-(l + «, T j exp^-- 



6(s - AEl) 



(23) 

The formulas for the primary-secondary and secondary- secondary contributions can 
be simplified to some extent but in the end some summations remain which must 
be done numerically. 

The number of nucleon-nucleon collisions as a function of s are plotted in Figure 
4. Both w = 0.5 and 0.8 are shown; there is little difference. The laboratory beam 
energy is 11.6 GeV per nucleon and the nuclei are gold. The spike represents the 
delta function from primary-primary collisions. The contribution from primary- 
secondary collisions falls from about 11 to 7 GeV -2 as s goes from 9 to 25 GeV. 
The contribution from secondary-secondary collisions is almost negligible. 

The pion multiplicity distribution arising from these hard collisions is shown in 
Figure 5. It drops by more than nine orders of magnitude in going from 6 pion 
production to 18 pion production. Typically there is only one hard nucleon-nucleon 
collision leading to the production of seven pions in a central gold-gold collision at 
this energy. 

We are interested in the possibility that one of these star bursts nucleates quark- 
gluon plasma. The precise criterion for this to happen is not known. However, 
we can make some reasonable estimates. In |4| we estimated that a critical size 
plasma droplet at these temperatures and baryon densities would have a mass of 
about 4 GeV. Any local fluctuation more massive than this would grow rapidly, 
converting the surrounding superheated hadronic matter to quark-gluon plasma. 
A similar estimate, based on the MIT bag model, a simpler hadronic equation of 
state (free pion gas), and with zero baryon density, was obtained much earlier ITT 



Another estimate is obtained by the argument that at these relatively modest beam 
energies most meson production occurs through the formation and decay of baryon 
resonances: A, N*, etc. The most massive observed resonances are in the range 
of 2 to 2.5 GeV. Putting two of these in close physical proximity leads to a mass 
of 4 to 5 GeV. We now need an estimate of the number of pions this critical mass 
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corresponds to. Let us assume that each particle, nucleon and meson, carries away 
a kinetic energy equal to one half its rest mass. If a particle would have too great 
a kinetic energy then it might escape from the nucleon-nucleon collision volume 
long before its neighbors and so would not be counted in the rest mass of the local 
fluctuation. Taking 4 GeV, dividing by 1.5, and subtracting twice the nucleon mass 
leaves about 6 pion rest masses. So our most optimistic estimate is that one needs 
a nucleon-nucleon collision which would have led to 6 pions if it had occurred in 
free space. One might be less optimistic and require the production of 8 or 10 pions 
instead. 

In Figure 6 we show the total number N > of nucleon-nucleon collisions which 
would lead to the production of at least n crit pions. We may view n crit as the min- 
imum number necessary to form a nucleation site or plasma seed. If n CT - lt = 6 is 
the relevant number then there are on average 7 such nucleon-nucleon collisions per 
central gold-gold collision. If 8 or 10 are the relevant multiplicities then there is only 
one such critical star burst every 1 or every 25 central gold-gold collisions, respec- 
tively. These numbers vary somewhat with w; the numbers quoted are averages. 
Conservatively, we may conclude that the probability of at least one plasma seed 
appearing via this mechanism is in the range of 1 to 100% per central gold-gold 
collision at the highest energy attainable at the AGS. These probabilities are about 
one to two orders of magnitude greater than those estimated in Q on the basis of 
thermal homogeneous nucleation theory. 



5 Consequences for the Multiplicity Distribution 

The results of the last section confirm the possibility of producing quark-gluon 
plasma droplets in rare events at AGS. Once formed the droplets grow rapidly due 
to the significant superheating of the hadronic matter. This process was explored 
in H] where it was found that the radii of such droplets can reach 3 — 5 fm. Since the 
phase transition is occurring so far out of equilibrium we would expect a significant 
increase in the entropy of the final state. This could be seen in the ratio of pions 



to baryons, for example, or in the ratio of deuterons to protons ||12|| . Along with 
the increased entropy should come a slowing down of the radial expansion due to a 
softening in the matter, that is, a reduction in pressure for the same energy density. 
Together, these would imply a larger source size and a longer lifetime as seen by 
hadron interferometry Jl3| . 



In this section we study one of the experimental ramifications in detail. Specifi- 
cally, we look at the charged particle multiplicity distributions and investigate under 
what conditions one might be able to detect the rare events from the structure of 
this distribution. 

In Figure 7 we plot the ratio of entropy to total baryon number S/B for the 
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hadronic and quark-gluon plasma phase for fixed beam energies. Fixed beam energy 
means that initially both the energy density and the baryon number density of the 
system is given which then determine the corresponding entropies via the equation 
of state. We use the equation of state discussed in |4j] for all further calculations. It 
is helpful to consider two extreme and opposite scenarios. Either the matter stays 
all the time in the hadronic phase, or the matter has been completely converted 
to quark-gluon plasma by the time to an d only hadronizes later. The difference of 
the entropies produced in these two scenarios is given by the difference of the two 
curves in Figure 7. It represents an upper limit on the additional number of pions 
produced. Since the temperature is comparable to or larger than the pion mass the 
excess entropy is proportional to the maximum number of excess pions 

AiV_ 1 AS 

3 = . 24 

B 3.6 B y ' 

The number of additional negatively charged pions per baryon AN_/B is linearly 
related to the entropy difference AS determined from Figure 7. The result is shown 
in Figure 8 for central Au + Au collisions. At beam energies of 11.6 GeV/A we 
produce 0.33 additional negatively charged pions per participating baryon. This is 
an upper limit, and in reality we would expect less. 

These additional mesons might be visible in the charged particle multiplicity 
distribution which would have the form 

P n = (1 - q) P h n ad (AM + qP^(N qg ) . (25) 

Here q is the probability of finding a central event in which plasma is formed, -P^ ad 
is the multiplicity distribution for purely hadronic events with mean iVhad, and P^ g 
is the multiplicity distribution for events in which a plasma was formed with mean 
N qg . 

Experimentally one would expect to see a bump in P n at larger values of n. 
A structure like that was found in charged particle multiplicity distributions in pp 
collisions at the CERN |14], |15[] and Fermilab fl6| , |3~7p colliders. For energies larger 
then 540 GeV a shoulder develops in the multiplicity distribution, becoming more 
pronounced as the beam energy increases. It is assumed that this structure is due 
to the onset of minijets. It is definitely an indication of new physics. 

In Figure 9 we plot the charged particle multiplicity distribution for pp collisions 
at y/s = 900 GeV from the UA5 collaboration flM. For energies less then 500 GeV 



it was found that the distribution could be well described by a negative binomial 
distribution of the form 



P n (n,k) 




n/k 



1 + (n/k) 



[1 + (n/k)) 



k ' 



(26) 
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The parameter k characterizes the width of the distribution. For k — > oo we recover a 
Poisson distribution, the distribution with the smallest width. One can see from the 
figure that at 900 GeV a single negative binomial (NBD) cannot describe the data 
anymore. A double negative binomial (DNBD) of the form discussed in eq. ( p5|) on 
the other hand describes it very well. The question remains to what extent a similar 
analysis might be able to reveal rare events of quark-gluon plasma production at 
AGS. 

A rough criteria for the observability of such structure in distributions of the 
form (BSI) is 

— J = P / !\ad = pqg , 27 x 

Here iVbin is the number of observed central Au + Au collisions for which the central 
multiplicity of the bin is A qg . The right-hand side of eq. (p7|) is the magnitude of the 
rare events to the overall multiplicity, while the left hand side gives the statistical 
resolution. The assumption here is that q is small, so that at A^ qg we can use 
P n ~ P^ ad for the left-hand side. 

To obtain a feeling for the shape and applicability of eqs. (p6|) and fl27D we plot in 
Figures 10 and 11 different negatively charged particle multiplicity distributions as 
might be expected for central Au + Au collisions at AGS with -Ebeam = H-6 GeV/A. 
From || we obtained the mean for purely hadronic events to be Ahad — 145. 
This is slightly larger than the value Ah a d = 131 ± 21 cited in |§ for 355 ± 7 
participating nucleons since we are assuming that all 2A nucleons are participating 
in the collision. The result depicted in Figure 8 for the upper limit on the additional 
number of negatively charged pions produced per participating baryon allows us to 
deduce an upper limit of A qg = 193 on the mean for the events with quark-gluon 
plasma production. In Figure 10 we plot the negatively charged particle multiplicity 
distribution defined in eq. (|25| ) for different values of the probability q. We use 
Poisson distributions for p had and P qg and take the upper limit for rare events 
A qg = 193 as the mean for P qg . A shoulder develops for small q and becomes more 
pronounced the larger q is. In Figure 11 we fix q = 0.1 and investigate the effect of 
different values of the mean A qg of the distribution for events with some quark-gluon 
plasma production. If this mean is close to the mean of purely hadronic events we 
will only find some broadening of the overall distribution. This would be the case if 
the phase transition is weakly first order or second order. For larger N qg we begin to 
see a well established shoulder develop. For large A qg a second maximum appears. 

It is clear that the exact values of the probability q and of the mean A qg of rare 
events will be crucial for the experimental observation of a phase transition. We 
have provided a first glimpse into this problem, but in the end it is up to experiment 
to discover new physics in multiplicity distributions at the AGS. 
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6 Summary and Conclusion 



We have estimated the probability that hard nucleon-nucleon collisions initiate the 
formation of seeds of quark-gluon plasma at AGS energies. Based on our previous 
studies we know that these will grow rapidly to convert most of the superheated 
hadronic matter to quark-gluon plasma. Our estimates are based on reasonable 
assumptions and approximations to the kinetic theory of hadronic physics. Better 
estimates could be made using event simulators like RQMD and ARC together with 
more detailed knowledge of multi-particle production in nucleon-nucleon collisions. 
We find that anywhere from 1% to 100% of central Au + Au collisions should lead 
to significant quark-gluon plasma formation. A major assumption is that there is a 
phase transition and that it is first order. 

We have already proposed that the formation of plasma in rare events should 
have an observable consequence for hadron interferometry, deuteron production, and 
the meson multiplicity distribution. In this paper we have studied the effect on the 
multiplicity distribution. It would be observable as a shoulder or second maximum 
at some multiplicity higher than the most probable one. If there is a phase transition 
but it is second order or weakly first order then the effect will be much more difficult 
to see. We eagerly await the results of experiments. 
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Figure Captions 



Figure 1: Schematic of a central collision between two nuclei. Nucleus 1 is incident 
from the left and nucleus 2 is incident from the right. The shaded area is the hot 
and dense overlap zone which is expanding along the beam axis with the original 
beam velocity v. A hard nucleon-nucleon collision leading to a large energy release, 
or star burst, is indicated at the longitudinal position z. 

Figure 2: The density of hadrons n tot in the hot and dense overlap zone as a function 
of laboratory beam energy E heara . 

Figure 3: The pion production cross sections cx n versus energy s as computed ac- 
cording to the text. Note that they are averaged over initial state isospin. 

Figure 4: Distribution dN/ ds in s of hard nucleon-nucleon collisions taking place in 
the hot zone. The value of w is 0.5 (4a) and 0.8 (4b). 

Figure 5: Number of hard nucleon-nucleon collisions N n leading to a particular final 
state pion multiplicity n. 

Figure 6: Number of hard nucleon-nucleon collisions N > with at least n cr - lt pions 
produced. 

Figure 7: Ratio of entropy to baryon number S/ B for fixed beam energy. 

Figure 8: Upper limit on the additional number of negative pions produced per par- 
ticipating baryon AN-/B in central Au + Au collisions as a function of beam energy. 

Figure 9: Charged particle multiplicity distribution for pp collisions at y/s = 900 



GeV. The parameters for the fits are taken from [15] with a probability of the sec- 
ond, high multiplicity, component being 0.35. 

Figure 10: Negatively charged particle multiplicity distribution for central Au+ Au 
collisions at i?beam = 11.6 GeV/A for different values of the probability q. The mean 
for purely hadronic events is taken to be A^ad = 145 while the mean for events with 
quark-gluon plasma production is taken to be N qg = 193. 
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Figure 11: Negatively charged particle multiplicity distribution for central Au + 
Au collisions at -Ebcam = H-6 GeV/A for different values of the mean N qg for rare 
events. The probability is fixed at q — 0.1 and the hadronic mean multiplicity is 
fixed at A" had = 145. 
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